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Abstract. We report a novel experimental technique that measures simultaneously 
in three dimensions the trajectories, the translation, and the rotation of finite size 
inertial particles together with the turbulent flow. The flow field is analyzed by 
tracking the temporal evolution of small fluorescent tracer particles. The inertial 
particles consist of a super-absorbent polymer that renders them index and density 
matched with water and thus invisible. The particles are marked by inserting at 
various locations tracer particles into the polymer. Translation and rotation, as well 
as the flow field around the particle are recovered dynamically from the analysis of 
the marker and tracer particle trajectories. We apply this technique to study the 
dynamics of inertial particles much larger in size (R p /r] ~ 100) than the Kolmogorov 
length scale 77 in a von Karman swirling water flow (R\ ~ 400). We show, using 
the mixed (particle/fluid) Eulerian second order velocity structure function, that the 
interaction zone between the particle and the flow develops in a spherical shell of width 
2R P around the particle of radius R p . This we interpret as an indication of a wake 
induced by the particle. This measurement technique has many additional advantages 
that will make it useful to address other problems such as particle collisions, dynamics 
of non-spherical solid objects, or even of wet granular matter. 
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1. Introduction 

Particle-laden flows are prevalent in natural and technological flows. For example they 
are relevant to warm cloud dynamics [U E] with impact on the climate, or to dust and 
pollution transport in atmosphere and oceans, or can be found in most technological 
processes where matter is mixed in fluids. In the past few years, significant advances 
have been made in experimental approaches to particle dynamics in turbulent flows 
thanks to the rapid development of visualization based three dimensional measurement 
techniques such as Lagrangian Particle Tracking (LPT) [HI [4]. When studying the fluid's 
dynamics with LPT, one has to assure that the particles behave passively and do not 
perturb the flow. This goal can be considered reached when the particles are density 
matched to the fluid, which avoids buoyant and inertial forces, and when their sizes 
are much smaller than the smallest length scale of the velocity gradients in the flow. 
Such particles are then called tracer particles. These requirements are very stringent 
and can be hard to achieve in some flows, e.g. thermal turbulent convection, where fluid 
density changes strongly in the thermal boundary layer 0, El [7j; highly turbulent flows, 
which exhibit very fast dynamics (fractions of milliseconds) on length scales that reach 
micrometers [8] (indeed L/rj oc i?e 3 / 4 , where L is the large scale of the flow, rj is the 
characteristic length-scale of the smallest eddy in the flow, also called the Kolmogorov 
length scale, and Re is the Reynolds number of the flow); or in the study of superfluid 
turbulence where the normal and superfluid coexist [9] . When the density of the parti- 
cles is different from that of the fluid, or the particle size is large, the particles cannot 
be regarded as tracers of the fluid and are called inertial particles. 

In this paper, we focus on the study of large particles in turbulent flows, i.e. 
particles much larger in size than the Kolmogorov length scale rj of the flow. The 
particles that we consider here have a radius R p ~ IOO77 ~ 0.1L. Experiments on 
two-particle statistics of small and heavy particles can be found in [TUl CD]- Most 
theoretical and numerical simulations of the complex coupling between the motion of 
inertial particles and their carrier fluid rely on simplified equations such as the Maxey- 
Riley-Gatignol equation p~2j,[l3]. This equation is only valid for point-like particles and 
couples the motion of the particle to the fluid, without any back-reaction of the particle 
motion on the flow (one-way coupling). This approach has been proven to be sufficient 
to describe the rich dynamics of small and heavy particles [H], but has also found its 
limit p~Ql HSl [HI (TTJ [HI IE] (even when corrected by the so called Faxen corrections) 
when considering particles larger than 477. To go beyond this simplified model, Hohmann 
and Bee [20j developed a direct numerical simulation using a dynamical pseudo- 
penalisation technique in order to satisfy the non-slip boundary conditions at the surface 
of a unique particle evolving in a turbulent flow. This strategy is very promising (as 
the one developed in [21], in which the particle is static) but also comes with high 
computational cost that currently restricts it to very low Reynolds numbers (R™ ax = 72 
for [20] and R™ ax = 20 for [21], where R\ is the Reynolds number based on the Taylor 
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(a) Side view (b) Top view 



Figure 1. Picture of the Von Karman water flow. The green sphere at the center of 
the apparatus represents the measurement volume. 

micro-scale of the flow). Therefore, for the foreseeable future, the investigations on large 
particle dynamics in highly turbulent flows must rely heavily on experimental work. Few 
experiments such as [22] have been conducted in 2D flows, also 2D studies of 3D flows 
based on Particle Imaging Velocimetry can be found in [23j , but a complete experimental 
3D approach resolving the particles dynamics and the flow has to our knowledge not been 
reported. The recent work of Zimmermann et al. [2H [25] has been able to resolve the 
full motion (translation and rotation) of single finite-sized particles in a highly turbulent 
flow. Here, we present a technique that goes beyond and resolves simultaneously the full 
three dimensional motion of the particles together with the turbulent flow field around 
it. 

2. Experimental setup 

In order to generate a highly turbulent flow in a relatively small experiment, we used a 
von Karman water flow, sometimes nick-named the "french washing machine" in the lit- 
erature (see Fig. CD) [8], [26], [27] . This flow is generated by two counter-rotating propellers 
of 28 cm in diameter facing each other. The propellers are driven at the same constant 
rotation rate. The turbulence chamber, shaped as a octagonal cylinder, measures 40 cm 
along the axis of the propellers and 38 cm in both height (vertical) and width (horizon- 
tal) in the cross-section. This apparatus is particularly well suited for LPT studies [8] 
since at the geometric center the mean flow velocity is zero. Therefore, at the center 
the particle dynamics is mainly driven by the turbulent velocity fluctuations (the large 
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Table 1. Parameters of the experiment, u' is the root-mean-square of the velocity, 
e is the energy dissipation rate per unit mass. L = u /3 /e is the integral length 
scale, r\ = (z^ 3 /e) 1 / 4 and = (yje) 1 ! 2 are the Kolmogorov length and time scales, 
respectively, where v is the kinematic viscosity of the fluid. Nf is the frame rate of the 
camera, in frames per and R\ = (lbu'L/v) 1 / 2 is the Taylor scale Reynolds number 
of the flow. 

scale flow being deterministic^) (2HI EZ] • 

To study the turbulent flow field, we seeded the water flow with red fluorescent 
polymer micro-spheres with a diameter of 107 pm and a density of p tr = 1.05 g/cm 3 ^. 
These particles behave as tracer particles since in our flow the Kolmogorov length scale 
IS fj 100 pm (see Tab. CD) and the fluid density is pf — 1 g/cm 3 . We measured 
three-dimensional particle trajectories with high spatial and temporal resolutions using 
LPT (3j [4] with three high-speed CMOS cameras (Phantom V10, manufactured by 
Vision Research Inc., Wayne, USA). The particle velocities and accelerations were 
then obtained by smoothing and differentiating the trajectories [28]. The measurement 
volume of the LPT was determined as the largest sphere that fits inside the complex 
volume defined by the intersection of the fields of view from all three cameras and the 
expanded laser beam. In this experiment, the sphere was 7 cm in diameter slightly 
smaller than the integral length scale L = (9.3 ± 0A)cm. The frame rate was set to 
2.9 kHz with a resolution of 768 x 768 (each pixel corresponds roughly to 100/xm in 
space). Given the massive amount of data required to conduct this experiment, we used 
the weighted averaging algorithm (simple and efficient) to locate the particles onto the 
2D images. The number of particles tracked per frame was of the order of 200. In 
this configuration, as demonstrated in [3], the 2D particle finding algorithm and 3D 
stereoscopic reconstruction processes ensure that the particle positions are recovered 
with an accuracy of 0.1 pixels « 10 pm. The relatively low seeding density further 
guaranteed that, using a "3 frames minimum acceleration" scheme [3], we obtained long 
and high quality particle tracks. 

We measured u r the velocity fluctuation and determined the energy dissipation 

X The ratio between the average velocity and the velocity fluctuations r = (wi)/v ((ui — (ui)) 2 ) does 
not exceed 25% over the entire measurement volume. 

§ These fluorescent particles are hard-dyed (internally-dyed) polymer particles which utilise the Firefli 
process to incorporate the dye throughout the polymer matrix. This method produces bright fluorescent 
colors, minimises photo-bleaching, and prevents dye leaching into aqueous media. They are made of 
polystyrene and are sold by Thermo scientific. Their absorption (emission) wavelength peaks around 
542 nm (612 nm). 
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(a) (b) 

Figure 2. Super- absorbent polymer particles, (a) Once immersed in water and in 
their dry form, (b) Beaker full of these particles crossed by a laser beam that is not 
refracted. The red luminescent dot in the center of the picture is a marker particle 
grafted at the surface of the solid particle (symbolized by the red dots in (a)). 

rate per unit mass e from the inertial range scaling of the second-order longitudinal 
and transverse Eulerian velocity structure functions from tracking tracer particles. As 
explained in details in [10], we also checked the consistency of the value using two exact 
inertial range relations: the Kolmogorov's "four- fifth law" and a theorem on the velocity- 
acceleration mixed structure function [29J, [30], EI] . The properties of the turbulent flow 
field are summarized in Table ED 

The aim of the investigation was to simultaneously measure the 3D trajectories of 
the tracer particles and the full dynamics and trajectories of finite size, inertial parti- 
cles. The most straightforward way was for the particles to be invisible in the fluid and 
to mark them with tracer particles. This avoided optical distortions and shadowing. 
The spherical particles were made of a super-absorbent polymer (poly-acrylate) from 
Aqualinos, which in their dry form had a diameter between 1 and 2 mm; once immersed 
in water their average diameter grew to 1 cm ~ L/10 ~ IOO77 (see Fig. [2] (a)). Thus the 
gel particles were 99.9% water and had almost the same refractive index and density 
as water. As shown in Fig. [2] (b) they are indeed almost invisible in water and do not 
refract the light path. 

In order to be able to observe the big particle's motion with the cameras, we grafted 
their surface with the same fluorescent particles that we used to track the turbulent flow 
field. As we will see later, we needed at least 4 of these markers per inertial particle in 
order to recover the particle's center and radius. Therefore, we injected them with 6 to 
8 markers, as shown by the red dots in Fig. [2] (ajjjl To achieve this we employed a lancet 
needle normally used by diabetics. We dipped it into the tracer particle powder to collect 
some of them on its tip and penetrated the surface of the big particle. After removing 
the lancet, tracer particles were left behind within the polymer gel. We quantitatively 

|| The regular separation presented on Fig. El (a) is not necessary, therefore we do not enforce it. 
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checked that the markers were indeed grafted to the big particles by injecting some 
of them exactly as depicted in figure [2] (a) thanks to a specifically designed injection 
guide. No relative nor absolute displacement could be detected from this well defined 
injection pattern. We marked about 150 big particles, which we inserted into the turbu- 
lent flow. This particle number was necessary to achieve that statistically, at any given 
time, at least one (and no more than two) big particles were in the measurement volume. 

The big solid gel-particles were elastically deformable, but the forces (shear and 
pressure) applied by the fluid onto them were not sufficient to deform them or trigger 
any internal flows (the internal effective viscosity of the gel being much greater than 
viscosity of the water). In fact, from measuring marker positions no deformation could 
be detected within the experimental uncertainties, while the particles were carried by 
the turbulent flow. Therefore, these gel-particle can, for the purpose of this study, be 
considered as solid. To prevent the particles from being damaged by the propellers, 
we placed grids (of 2x2 mm 2 mesh size and 75% aperture) in front of both of them 
(see Fig. []]). Measurements showed that the influence on the flow was negligible. We 
ensured that the measurement volume was sufficiently far from the grid so that the 
large particles' dynamics was not influenced by the meshes. This can be seen from the 
following analysis: the viscous relaxation time of the particle is, 



where p p and pf are respectively the particle and fluid density, d p the particle diameter 
and v the kinematic viscosity of the fluid. This applies only when the particle Reynolds 
number (Re p = u'd p /v) is very small. Following J32], [33] , one can show that the relevant 
time scale that takes into account small but finite particle Reynolds number can be 
written as: 



where the exponent n also depends on Re p : n = 0.82 — 0.051og 10 Re p . With p p /pf ~ 1, 
d p = 0.01 m, v = 10 -6 m/s 2 , and v! = 0.1 m/s this time scale is r p = 0.38 s, which yields 
L p = u'r p = 3.8 cm beyond which the big, inertial particle has reached the fluid veloc- 
ity. The minimal distance between the grid and the measurement volume is ^ 7 cm on 
both sides and thus we expect that the particles' dynamic is not impacted by the screens. 

As pointed out before, we were seeding the turbulent flow with two types of 
particles: big particles marked by fluorescent tracers and the fluorescent tracers 
themselves. When tracking all fluorescent particle tracks where visible simultaneously; 
marker particles trapped in the big particles did not change their separation, while 
tracer particles would separate quickly. 




(i) 
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3. Identifying the finite-size particles 

One very significant advantage of the novel technique presented here is that it allows to 
follow two species of particles using only one LPT system. For comparison, Guala et 
al. [34J used two LPT systems, which makes such an experiment very demanding. To 
identify the two populations of particles we took advantage of the fact that the flow was 
highly turbulent. It is well known that tracer particles separate very quickly [35], [TO]. In 
the contrary markers fixed to the surface of a big particle will stay at a constant distance. 
Specifically, two tracer particles initially separated by a distance R = |R(t = 0)| in the 
inertial range (i.e. rj <C Ro <C L) will separate on average following a ballistic regime 
predicted by Batchelor [36] such as: 



for t <C t = (Rl/e) 1 / 3 , where SH(t) = R(t) — R(t = 0) is the vectorial separation 
increment and C 2 = 2.1 as suggested from a compilation of available data [37] . 

3.1. Finding pairs of tracks 

The first step is to find pairs of fluorescent particles tracks whose distance to each 
other stays constant during the time of observation. As the fluorescent markers are 
embedded in the particles with radius R p) we scan only for tracks whose separations 
are 0.05 < Rt/R p < 2(here R p = 8 mm). The lower bound comes from the fluorescent 
particles injection process that we used. Additionally, to limit possible false choices, 
we only conduct this analysis on pairs of tracks that coexist for longer than 2r v . 
This is the characteristic time-scale after which the acceleration components of a 
fluid particle decorrelated entirely [8]. To estimate the tolerance on the separation 
increments, we have used Equ. [3] with R = 0.05R P , t = 2r v and e = 0.011 m 2 /s 3 which 
gave us \J (5H 2 ) w 0.1R P . In the analysis presented here we used a slightly tighter 
criterion: \J (5H 2 ) < 0.01i? p , that corresponds approximately to our spatial resolution 
(100 jim/ pixels). 

3.2. Finding groups 

To ease the explanations, two tracks that were paired by the algorithm described above 
(section 13.11) are now called " friends" . We can now identify the fluorescent particle 
tracks that belong to a unique inertial particle. To do so we simply apply the adage 
"the friend of a friend is my friend". This allows to generate groups of tracks that trace 
the inertial particles' motion. By additionally imposing that at least four tracks coexist 
for more than 2r v we exclude optical artifacts. Such a group is presented on Fig. [3] (b). 




(3) 




Figure 3. Track reconnection inside the groups, (a) Definition of the search volume 
for possible track reconnection. (b) Group before reconnection. (c) Group after 
reconnection. (d) Probability density function of the track length within the groups, 
before and after reconnection. 



3.3. Track reconnection inside the groups 

As one can see from Fig. [3] (b), the tracks that belong to a group occur in short 
segments. This can be attributed to limitations of the experiment, which include 
variations of the illumination intensity, particles blocking the line of sight of the cameras, 
particles blocking the illumination beam (shadowing), the presence of light-insensitive 
circuits on the CMOS sensors of the cameras, and the effect of thermal, electronic and 
environmental noise (see [4]). To reconnect the tracks we use an algorithm similar 
to that developed by Haitao Xu at the Max Planck Institute for Dynamics and Self- 
Organization j4]. Within a group, at the end of each track, we define a cylindrical 
search volume, parametrized by its length Cx, its radius Cr and its axis oriented along 
Vi (the velocity at the end of the track-segment). In Fig. [3] (a) such a cylinder is shown 
exemplarily in green. Tracks of the group that have their starting point within this search 
volume are considered potential re-connection candidates (track 2 on Fig. [3] (a)). Then 
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for each candidate, we compute the linear interpolation of the trajectory and compare 
the beginning and ending interpolated velocities (resp. V 3 and V 4 on Fig. [3] (a)) to the 
actual measured ones (resp. Vi and V 2 on Fig. [3] (a)). If the beginning and ending 
relative velocity differences are smaller than a tolerance Ct (e.g. |V 3 — Vi|/|Vi| < Ct) 
we consider these as a possible reconnection. If there are several candidates, we re- 
connect the best match. We remind the reader that we are working within a group 
obtained by the previous step (section 13.21) . that is to say a very small amount of 
tracks compared to the work in [4J. In our experiments, we found that Cr = 0.5 mm, 
Cl = 30 mm and Ct = 0.5 give satisfactory results. For example, we present the result 
of this reconnection process on a particular group in Fig. [3] (c), and its statistical effect 
on several measurements in Fig. [3] (d) (please note that this more than doubles average 
track length inside the groups). 

3.4- Cleaning the groups 

The process described above is not error-proof. Indeed, while scanning examples of 
groups by eye we could identify wrong tracks. We therefore developed the "group 
cleaning" technique. In this step, we scan each group and we keep only tracks that 
have more than three direct "friends" as defined in section 13.11 We use the same 
parameters except that we relax the constraint on the maximal separation increment 
to ^(STL 2 ) < 0.07R P . This we can do as the data was already connected by linear 
interpolation. 

3. 5. Concluding remarks on the particle tracking procedure 

As shown by Fig. [4] we can now separate the two populations of fluorescent particles. 
In this process, as discussed above, we interpolate trajectories (section 13.31) solely to 
connect tracks. These interpolated data points need to be excluded when statistics of 
the fluid velocity, or vorticity along tracks is considered. We would also like to note that 
an inspection by eye of « 5000 groups did show the procedure described above works 
only well, when only one inertial particle is in the measurement volume at any given 
time. We therefore discarded manually the rare cases where two or more particles (some 
colliding) were in view. This did not limit the statistics significantly. In the future it 
will be not hard to modify the algorithm to include this automatically. 

4. Translation of the finite size particle 

Here, we describe how we extract the center position x c and radius R p of the big inertial 
particles for each group identified by the procedure described in section [31 These groups 
have a minimum of four marker-trajectories. As described above, our injection technique 
ensured reproducible and well defined penetration depth into the inertial particle. The 
injection was such that the markers were inserted not more than 0.5 mm from the 
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Figure 4. Result of the process described in section [3] on a 3 s movie, (a) All 
the trajectories of the fluorescent particles; (b) Same data set in light blue. The 
red trajectories correspond to the fluorescent particles bounded to the surface of the 
super-absorbent polymer spheres identified by the procedure described in section [3l 



polymer surface. The center position x c and radius R p of the inertial particle is entirely 
defined by the sphere equation: 



(xi - x c ) 2 + (yi - y c ) 2 + (zi - z c f 



Rl 



(4) 



where [x^y^Zi) are the z-th marker's 3D-coordinates. For each permutation of four of 
these markers, one gets the following linear system: 

P-{A B C D) T = Tr(X • X T ) 



(5) 



with 



p = 



V 



-1 xt y 1 zt \ 

-1 x 2 y 2 z 2 
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(6) 



and A = x 2 c + y 2 c +z 2 c 



( x\ 2/i z\ \ 

X 2 V2 Z 2 

%3 VS z?> 

\ V4 Z4 ) 

R 2 . B = —2x C) C = — 2y c and D = — 2z c . Let us consider a group 
that contains n > 4 markers trajectories at a given time t, the number of permutations 
of 4 trajectories out of n is M = ^0z^j\ = (4)- Therefore, at each time step, we solve 
M times the linear system given by Equ. [51 We define x C5m (i) the center coordinates 
obtained by the ra-th permutation and d m (t) = Ylt=i !>2j=i+i l x i,m(i) _ Xj >m (£)| the sum 
of the distances between the markers. To define the actual center coordinates of the 
inertial particles, we perform an average of x C5m (i) weighted by d m (t)\ 

M M 

Xc(t) = ^ x c , m (t)d m (t)/ ^ d m (t). (7) 



m=l 



m=l 
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We have made this choice as to correct for a non-homogeneous distribution of the 
markers on the inertial particle's surface. Indeed, four markers homogeneously 
distributed around the sphere define a center position and a radius more accurately 
than four agglomerated points. 
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Figure 5. Results of the center finding process for one example group, (a) Trajectories 
(thin lines) identified to belong to the particle and the center (thick green line), (b) 
Radius (bold black line) and the distances from the center to each trajectory (solid 
lines) for the tracking time of the particle. In this specific case, the particle was tracked 
for 1115 frames = 41.3^ « 0.4s). The inset at the bottom shows the 3 coordinates of 
particle center x c subtracted by its filtered version (Gaussian kernel of a = 3.5 frames) 
xj, in [mm]. 



The result of the center finding process on a particular group is shown on Fig. [5] (a), 
together with the radius of the inertial particle and the distances from the center 

position to the markers of this group r^(i) = |xi(i)— x c (i)| in Fig. 0(b). The latter shows 
that we are able to measure the particle radius with an uncertainty of ±200/xm, which 
corresponds to the uncertainty of the fluorescent marker insertion. We also observed 
that R(t) exhibits unrealistic rapid variations when a marker track appear or disappear 
from the group (around frames 200, 350, 800 and 900). These kinks are non-physical 
and result from the uncertainty of the marker position inside the inertial particle. 
Just as for tracer particles [28], we measured the velocity of inertial particles from the 
center trajectory by convolving with a properly normalized, truncated, differentiating 
Gaussian smoothing kernel. The characteristic time scale of the filter was chosen as 
a = 3.5 frames = 0.13^ sufficient to minimize the effect of noise on the tracer particles' 
acceleration variance. At the kinks observed in the center position of the big particles 
x c (t) (Fig. [5] (b)) the values of the velocities were very large, discontinuous, and thus 
unphysical. We flagged and excluded those from the analysis together with 2a = 7 
frames before and after those discontinuities. 
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5. Rotation of the inert ial particle 

To measure the rotation of the inertial particles, we used the center position together 
with the trajectories of the markers around the center of the sphere. To do so, we took 
advantage of an algorithm first introduced by W. Kabsch [38j [39] that was developed 
to compare molecular conformations in chemistry. This algorithm seeks the optimal 
rotation matrix U between two sets of points by minimizing the root mean square of 
their separation. 

Let us consider two sets x and y of N paired points that have the same centroid at the 
origin. We are looking for the 3x3 unitary matrix U that aligns best x with y. This 
can be achieved by minimizing the root mean squared: 

1/2 

(8) 

That can be rewritten as: 

N 

ND 2 = J2[( Ux >) 2 + Vi ~ 2Ux iVil ( 9 ) 

1=1 

where one can see that minimizing the left hand side is equivalent to maximize the 
negative term on the right hand side. If we represent the points of set x (resp. y) by a 
TV x 3 matrix X (resp. F), the quantity to maximize is: 

N 

Ux iVi = Ty(Y t UX) = Tt((XY t )U) (10) 

t=i 

where XY T is a square 3x3 matrix, that can be rewritten using its singular value 
decomposition (SVD) as XY T = V SW T . V and W T are orthogonal matrices of the 
left and right eigenvectors of XY T and S is a diagonal 3x3 matrix containing the 
eigenvalues s\ > s 2 > 53. Using the commutation properties of the trace operator, 
Equ. [10] can be rewritten as: 

3 

Tt(VSW t U) = Tt(SW t UV) = ^ SiTa (11) 

i=i 

where T = W T UV is a product of orthogonal matrices and itself an orthogonal matrix 
with det(T) = ±1. Therefore each elements of T are equal or smaller than 1. Thus to 
maximize Eq. [IH T has to be the identity matrix T — I. 
The last constraint that we have to consider is that 

U = WTV T (12) 

has to be a proper rotation matrix with det(U) = 1 (rows/columns of U have to form 
a right handed system) [40]. For det(Xy T ) > 0, it follows that det(U) = 1, but when 
det(XY T ) < 0, det(U) = —1. In the later case, one has to settle for the second largest 



D 



1 N 
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value of Eq. [TTJ This value is obtained when T 33 = — 1 since si > s 2 > s 3 . This finally 
allows us to write the optimal rotation matrix U as 





(l 










u = w 





1 














d 


) 


where d = sign(det(XY T )). 









When applying this algorithm to our data, we constructed the N x 3 matrices 
X n i = x ni (t) — x Cj i(t) and Y ni = x ni [t + At) — x c ^(t + At), where x ni are the position 
components of the n-th marker trajectory that existed at times t and t + At (At is 
equal to the acquisition time) . By construction for each group N > 4 (see section Ej) 
thus ensuring a proper rotation matrix. Then for each time steps, we computed the 
covariance matrix C = XY T , its SVD C = VSW T and the sign of its determinant 
d = sign(det(C)). Next, we applied Equ. [13] to obtain the optimal rotation matrix U. 




Figure 6. (a) Reduced probability density function of the rotation rate component of 
the inertial particle (the different symbols corresponds to the 3 components uo Xl uo yi uj z ). 
The dashed line corresponds to a stretched exponential function, (b) 3D visualization 
of the inertial particle trajectory together with the tracer particles around it. The 
colors encode the velocities of the particles, warm colors indicate fast particles and 
cold ones low velocity particles. 

From the rotation matrix, one can extract the rotation vector or the Euler angles 
and compute the rotation rate of the inertial particle cj. One can estimate an uncer- 
tainty of about 5% on the components of uo from the analysis of the residuals of the 
least square optimisation described above. Figure [6] (a) shows the measured probabil- 
ity density function of the components for the rotation of the inertial particle. This 
distribution is non Gaussian and can be described by a stretched exponential function 

with a best fit parameter s = 0.55. Zimmer- 
form with s = 0.45 which is in good agreement 



p(„\ - exp(3s 2 /2) 

mann et al [21] used the same functional 



1 _ erf / \n(\x/V3\)+2s 2 
\/2s 
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with our results. 

As shown in Fig. [6] (b), we measured simultaneously the full motion of a finite size 
particle, translation and rotation, and the highly turbulent flow field from tracer particle 
tracking. This method thus allows us to study experimentally the interactions between 
the particle and the flow. 

6. Particle flow correlation 

This novel measurement technique is particularly well suited to study particle fluid 
correlation. The quantity that we want to study here is the second moment of the 
longitudinal velocity difference between the inertial particle and the flow, also called the 
mixed (particle/flow) Eulerian second order longitudinal velocity structure function: 

^r(r) = ([(v P (0)-v f (r))-f] 2 ) (14) 

where v p and Vf are the inertial particle and fluid velocity vectors, and r is the separation 
vector between the center of the inertial particle and the tracer particle considered 
(f = r/|r|). In general, this quantity depends on the position of the particle and on the 
separation vector r. Here however, we measured at the center of a von Karman mixer, 
where the flow is fairly homogeneous and isotropic, therefore one can assume that D™[ x 
is only a function of r. Additionally, given the symmetries of the system, we averaged 
the data azimuthally around v p axis, where we expect rotational invariance. 

In order to quantify the influence of the inertial particle on the turbulent flow field, 
we have normalized D™[ x by its expected value D^™ in the case of a point-like particle 
(a tracer/fluid particle since p p /pf ~ 1): 

DTl = nTf/DO? = Dr L x /C 2 (erf/ 3 (15) 

where D^™ = (^(er) 2 / 3 is the second order longitudinal velocity structure function of 
a homogeneous and isotropic turbulent flow (neglecting intermittency corrections) for 
separations r in the inertial range (77 <C r <C L). 

D™ x is shown in Fig. The statistical convergence of this quantity, obtained by 
binning the space around the big particle, was checked by splitting the dataset in two 
sub-samples, leading to the same conclusions detailed hereafter. First, one can see that 
both limits r = R p and r — > 00 are consistent with what we would expect. Indeed, 
D™[ x (r = R p ) = simply means that no fluid is entering or exiting the solid particle, 
and D™l x (r — > 00) = 1 shows that the particle/flow interaction occurs locally. We 
can identify a spherical shell of radius 1 < r/R p < 3, where the correlation between 
the inertial particle and the flow is greatly enhanced. The numerical results by Naso et 
al. [21] are in good agreement with our experimental observations. Furthermore, beyond 
an isotropic spherical shell of distance R p from the surface of the inertial particle we 
observed a weak anisotropy with respect to the particle direction. This could indicate 
the presence of a wake induced by the particle. 
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Figure 7. Normalized mixed longitudinal 2nd-order Eulerian velocity structure 
function D™L X • ^'l'l as a function of the relative separation distance r/R p (where R p 
is the particle's radius) and being the angle between the separation vector r and the 
velocity vector of the ball v p . The orientation of v p is shown by the grey arrow (with 
an arbitrary length). The grey circle in the middle of the map symbolizes the inertial 
particle. Assuming axis- symmetry around the velocity axis of the particle, the lower 
half is a reproduction of the upper half mirrored on the center line. The black dashed 
lines show the separations r = 2R pi 3R P , 4:R pi 5R P . 

7. Conclusion and perspectives 

The main advantages of this new measurement technique are: 

(i) it requires nothing more than the usual equipment used for "standard" LPT or 
Tomo-PIV (3 or 4 fast cameras and a light source); 

(ii) it allows to follow in three dimensions the full dynamics (translation and rotation) 
of a solid body (here a spherical particle) together with the flow field carrying it; 

(iii) the solid object, being made out of super-absorbent polymer, has the same index 
of refraction as the fluid which makes it invisible and therefore it does not block 
the fields of view of the cameras (no shadowing effects); 

(iv) the geometry of the solid object is recovered during the acquisition (in this case the 
particle radius R p ) and could be studied dynamically; 

(v) several solid objects can be followed simultaneously in the measurement volume. 

This list demonstrates the capabilities of this novel measurement technique. Here 
we applied it to the dynamics of big inertial particles in turbulent flows (see section 6) 
and identified the zone of interaction between the particles and the flow. Other questions 
that can be addressed are, for example, the study of collisions of particles with equal and 
different particle sizes, and the dynamics of non-spherical objects in turbulent flows. At 
high particle number density this technique can also be used to investigate the dynamics 
of wet granular matter and of the flow between the grains. 
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